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Abstract 

We investigate the consequences of the pseudo-complex General 
Relativity within a pseudo-complexified Roberston- Walker metric. A 
contribution to the energy-momentum tensor arises, which corresponds 
to a dark energy and may change with the radius of the universe, i.e., 
time. Only when the Hubble function H does not change in time, the 
solution is consistent with a constant A. 

PACS: 02.40.ky, 98.80.-k 

1 Introduction 

In the past several paths have been proposed on how to extend the 
theory of General Relativity (GR). One of the first was Einstein himself 
[H [2] who, in an attempt to unify electrodynamics with GR, extended 
GR to complex GR. For more recent articles on complex GR, you may 
consult [31 m . Others [3 El [3 El El [TOl HH H21 H31 [IH [J5] introduced 
a maximal acceleration, related to a minimal length parameter in the 
theory, or [161 El HS1 H] tried to extend GR by introducing hyperbolic 
coordinates. 
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The introduction of a maximal acceleration and the use of hy- 
percomplex (a synonym for para-complex) coordinates is very much 
related to what was published in [20]. In [20] a pseudo-complex (pc) 
extension of General Relativity (GR) was proposed. (To the article 
[20] we will refer from here on as (I).) The main objective was to inves- 
tigate the analogue of the Schwarzschild metric and its consequences 
within the pseudo-complex extension of GR. Its possible experimental 
verification was and is of prime importance. The pc-GR is formulated 
in a different manner than in [T7\ FL8\ [T9] . Because pseudo-complex 
numbers exhibit a zero divisor basis (see section II), we were able to 
define two independent theories of GR in each zero-divisor component, 
which were later connected. The resulting length squared element is 
similar but not equal to |3 El El El El HDl IIH US H3 [131 H5I ■ It contains 
an additional term proportional to the velocity times acceleration. As 
already said, the analogue of the Schwarzschild solution was studied 
within the new formalism. A strong heuristic assumption was made, 
requiring that the pc scalar curvature 7?- = 0. This did lead to strong 
corrections, which are excluded by the Parametrized-Post-Newton for- 
malism (PPN) [21] . However, relaxing the condition gives corrections 
which are smaller and not contradicting the PPN formalism; details 
will be given in a forthcoming publication. Deviations of the redshift 
to the standard Schwarzschild case were calculated. 

One of the main results and messages in (I) was: The pseudo- 
complex description, with its extended variational principle, contains 
no singularity, i.e. black holes don't exist. Furthermore, the new the- 
ory (I) automatically introduces dark energy, whose density depends 
in the analogue of the Schwarzschild case on the radial distance. In 
fact, the increase of the dark energy density towards smaller radial 
distances finally hinders (prevents) the collapse of a large mass. The 
increase of dark energy, as a function of the radial distance, depends 
sensitively on some functions which are very difficult to deduce. 
The results are in line with the model described in [22] where the dark 
energy is described by a scalar field coupled to gravity. It is interesting 
and quite helpful for understanding, to compare our results with those 
obtained within this model. In the model |22] the relativistic equa- 
tions are solved numerically. As a result, dark energy accumulates 
around the central mass, reducing the radius of the event horizon. 
The advantage of this model [22] is that it provides a distribution of 
the scalar field, whose intensity should be proportional to our £ func- 
tions. In one case, the dark energy falls off very quickly with the radial 
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distance. The disadvantage is the large numerical effort needed. The 
similarity between this procedure and our theory, however, indicates 
that both investigations might be useful for each other. 

Here, in this contribution, we intend to show that the same effect 
introduces a dark energy in models of the Robertson- Walker (RW) 
type universes [29\ 128] , Our theory also contains a minimal length 
scale. Its influence on solutions might be important for large mass 
concentrations [23]. Here, we will not discuss it but refer to a later 
publication. Whether our theory is realized in nature also depends on 
the predictions which can be made and on their experimental verifica- 
tion. In the present contribution our principal objective is to identify 
some consequences of the pseudo-complex theory, and to determine 
whether possible differences to Einsteins General Relativity may be 
observed. 

In section II a short review on how to define pc variables and their 
properties will be given. This is for those readers who are not yet 
accustomed to this mathematical structure. A simple presentation can 
also be found in (I). In section III we will formulate the pc extension 
of the Robertson- Walker metric. It will be quite straight forward and 
the steps can be copied from [29] or any other book on GR may be 
consulted. For details, we will give mainly reference to [29] . In section 
IV we will solve the pc-HW model and in section V we shall discuss 
some consequences. Section VI summarizes the main conclusions. 



2 Pseudo-complex variables and pseudo- 
complex metric 



Here we give a brief resume on pseudo-complex variables, helpful to 
understand various steps presented in this contribution. The formulas, 
presented here, can be used without going into the details. A more 
profound introduction to pseudo-complex variables is given in [244125]. 
which may be consulted for better understanding. 

The pseudo-complex variables are also known as hyperbolic \16\ [T7] , 
hypercomplex [26] or para-complex [27] . We will continue to use the 
term pseudo-complex. 

The pseudo-complex variables are defined via 



X = x\ + Ix 2 
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with I 2 = 1. This is similar to the common complex notation ex- 
cept for the different behavior of /. An alternative presentation is to 
introduce the operators 

(2) 

with 

o± = a± , cr + (T_ = . (3) 

The a± form a so called zero divisor basis, with the zero divisor defined 
in mathematical terms by P° = P°_UP°,, with P° ± = {X = Xa±\X e R}. 
The zero divisor generates all the differences to the complex number. 
This basis is used to rewrite the pseudo-complex variables as 

X = X + a + +X^a_ , (4) 

with 

X± = Xi ± X 2 

or 

xi = \(X + + X_) , x 2 = i(X+-X.) . (5) 
The pseudo-complex conjugate of a pseudo-complex variable is 

X* = xi- Ix 2 = X+a- + X-a + . (6) 
The norm square of a pseudo-complex variable is given by 

\X\ 2 = XX* = x\ - x\ = X+X_ . (7) 

This allows for the appearance of a positive, negative and null norm. 
Variables with a zero norm are members of the zero-divisor, i.e., they 
are either proportional to a + or cr_. 

It is very useful to carry out all calculations within the zero divisor 
basis, a±. Here, all manipulations can be realized independently in 
both sectors, because a+a- = 0. 
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In each zero divisor component, differentiation and multiplication 
can be manipulated in the same way as with real or complex variables. 
For example, we have [25] 

F(X) = F(X + )a + + F(X_)a- (8) 
and a product of two functions F(X) and G(X) satisfies 



F(X)G(X) 
= (F(X + )a + + F(X_)a_) (G(X + )a+ + 

= F(X + )G(X+)a+ + F(X-)G(X-)a- , (9) 

because <7 + cr_ = and o\. = a±. As a further example, we have 

F(X) F(X+) F(X.) 

G{X) ~ ~G{X+) ~G{xZ) ■ (10) 

This can be proved as follows: 



F(X) _ F(X + )o- + +F(X_)o-_ 
GjX) ~ G(X + )o- + +G(X_)cr_ 

{F{X + )a + +F(X.)a^)(G(X + )a^+G{X_)a + ) 
(G(X+) CT+ +G(X_) CT _)(G(X + ) CT _+G(X_) CT+ ) 



(11) 



where we have multiplied the numerator and denominator by the 
pseudo-complex conjugate of G(X), using a\ = With <r + cr_ = 
and o± = a±, the last expression can be written as 



(F(X + )G(X-)v + +F{X-)G(X + )a-) 
G(X + )G(X_)(a ++( x_) 



(12) 



Because a + + cr_ = 1, we arrive at Eq. (|T0 
Differentiation is defined as 



DFjX) F(X + AX)-F(X) 

~~DX~ = afeo AX ' (13) 

where D refers from here on to the pseudo-complex infinitesimal dif- 
ferential. 
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A very important difference to the standard GR is the introduction 
of a modified variational principle. It states that the variation of an 
action has to be within the zero-divisor, i.e. 

JSeP , (14) 

where "P denotes the zero divisor given by all values which are either 
proportional to <r + (A<7+) or to cr_ (A<7_). For convenience, the latter 
is chosen in this contribution, as it was in (I). However, instead of 
"e P " we write here "= A<r_". The number at the right hand side of 
(|14|) has zero norm and can be treated as a "generalized zero", thus 
representing a minimal extension of the variational principle. 

In the pc-GR the metric is pseudo-complex, i.e., (see (I) for details) 

9tiv = 9^ v cr + + g~ u a- . (15) 

Each component (o~±) can be treated independently, for many pur- 
poses, e.g., as how to define parallel displacement, Christoffel symbols, 
etc. Most of the steps, known from standard GR can be carried out 
analogously. This is the advantage of using the zero-divisor basis. It 
can also be shown that the four-divergence of the metric tensor is zero 
(see (I)). 

The connection of both components of the zero-divisor basis hap- 
pens through the variation principle, mentioned above. It states that, 
with the convention used and 5S = 55+cr+ + 5S-0~-, we have 

5S+ = and SS- = A . (16) 
For the length element square 

du 2 = g^DX^DX" 

= g+ v DX%DX\a+ + g~ v DXtDX"a + 

= DX+DXla + + DX-DX^a_ , (17) 

reality is imposed through 

du* 2 = du 2 . (18) 
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This implies that du 2 is equal to its pseudo real part 



du? = -(dx+dx1 + dx;dx^) 

= dx^dx^ + l 2 du ll du )1 (19) 

where in the last step we simply substituted X^ and X± by their ex- 
pressions in terms of the coordinates and four-velocities. The pseudo- 
imaginary component has to vanish, i.e., 

= ^DX+DX^-DX'DXt) 

= lidx^du 1 * + du fl dx^) . (20) 

One way to proceed is to take (|19p as the final length element and 
impose (120p as a constriction, noting that (120p gives the dispersion 
relation: Integrating, this gives UpU^ = 1. This is a nice feature 
because in other models the dispersion relation is usually put in by 
hand. 

Using Eq. (37) of (I), which reads in differential form 
dXfi = 9%dx v + lh^ u du u 

Idu^ = lg® u du u + h iiv dx v , (21) 

the length element (I19p and the constriction (I20p can be respectively 
written as 



du? = g^dx^dx" + l 2 du^du v ) 

+lh flu (dx' l du u + dvPdx v ) , 

(22) 



and 



hfu, {dx^dx v + l 2 du^du v ) 
+lg° v (dx^du u + du^dx v ) = . (23) 

These are the expressions reported in (I). 
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Some models [3j H], using complex coordinates, start with a real 
length element squared, defining it as \ [doj* + doj]. If one does the 
same with pseudo-complex coordinates, one misses then the simplifi- 
cations involved, treating first the zero divisor components separately. 
In the complex version [U 0] there is no zero divisor and, therefore, 
one does not have this advantage. 



3 Pseudo-complex Robertson- Walker 
Metric 



Previously we followed the steps in chapter 12.3 of the book of Adler- 
Bazin-Schiffer [29]. Of course any other book on GR can be consulted. 
In this section we shall mainly repeat these steps, for sake of complete- 
ness, with the difference that the variables are now pseudo-complex. 
The reader will see that the formulation is identical to standard GR, 
with the difference of the appearance of additional functions due to 
the modified variational principle. 

In order to proceed, one choses so called Gaussian coordinates, 
in which one uses a distinguished (absolute) time coordinate, thus 
the abandonment of a completely covariant treatment of the cosmo- 
logical problem [29]. This is the price one has to pay to simplify 
the cosmological models and to describe physical reality in convenient 
mathematical terms. Such coordinates were first introduced by Gauss 
within a different context. 

The pseudo-complex length element in Gaussian coordinates, be- 
fore imposing reality, is given by 



where we already used the pseudo-complex coordinates. G is a func- 
tion of time and the radial coordinate R. As will now be shown, 
G can be written as the sum of the functions g(X°) and f(R), the 
first one depending only on time and the second one only on R, i.e., 
G(X°,R) = g(X°) + f(R). The starting point is the equivalence 
principle that two observers at two different points observe the same 
physics. The only difference may be in the scale the two observers use. 



dco 2 = (dXy-e G{ * x0 ^(dR 2 + R 2 d0 2 + R 2 sin 2 ed<i) 2 ) 




(24) 



8 



Thus the ratio of the proper distance element at two different space 
points Ri and R 2 must remain fixed in time: 



const in time 



(25) 



e G(X°,R 2 ) 

i.e., this ratio must be independent of X°. Therefore one must have 



G(X°, Ri) = G(X°,R 2 ) + F(R 1 ,R 2 ) 
which then yields for Eq. (|25p 



(26) 



e G(X ,_R2)+F(Ri,i?2) 
e G(X°,R 2 ) 



3 F(_Ri,R 2 ) 



which is independent on time. 

If we choose a fixed value for R 2 we can write 



(27) 



Christoffel symbols: 

The equation for the geodesies is given by 



5 J [{X ) 2 - e G (k 2 + R 2 6 2 + R 



(28) 



2 sin 2 66 2 



ds e P L 



(29) 



with s being a curve parameter. In addition, we used the new varia- 
tional procedure, requiring that the variation gives a number within 
the zero-divisor basis. 

After variation, the following equations of motion are obtained (a 
dot refers to the derivation with respect to s, the curve parameter, 
and a prime indicates for the function g a derivative with respect to 
X° while for / it is a derivative with respect to R). 
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X° + \g'e G ( R 2 + R 2 9 2 + R 2 sinH 



R+\f'R 2 +g'X°R 

_ (1/' + ^ (^2 + R 2 sin 2 e ^2 

+ 2 (i/' + i?0 + g 'X°0 - sin6cos6<j) 2 
4> + 2 Uf + ±)R4) + g'X°<j) + 25000*0 



(30) 



We used the convention that on the right hand side stands an element 
in the zero divisor basis proportional to a_. Choosing it proportional 
to o" + would give equivalent results, i.e., the <r_ and er + components 
are just interchanged. 

Comparing this with the equation of motion 



X» + 



A* 

i/ A 



(31) 



yields the non-zero Christoffel symbols (others can be deduced using 
the symmetry properties of the Christoffel symbols): 





1 1 





1 

1 

1 1 
1 



2 y 



1 



~-g'e G R 2 



-g'e G R 2 sm 2 t 
T 



V 

2 



(]/' + 1 ) 
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2 

3 3 



= — sin#cos# 

= cote . (32) 
^From the line element we find for the determinant of the metric tensor 



3 

2 3 



lny^ = ^g(X ) + p(R) + 2lnR + ln\sm9\ . (33) 



Using the Christoffel symbols given in (I32p . one finds 



A* 


/' 



/' 



/' 



2 J ^ R J ^ R 2 



R 2 



1 



G il , J2 



e 5 +5 



2 J ^ 



-cos 2 



R 2 8in 2 e 



The following relations are also useful 



(34) 



i/ | I /i 



-/2 



1 i/ 



1 „ 
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/' 



/' 



-e G g' 2 



o/2 



R 



r 



R 2 



l e G n _l f f2_l f , 
2 y 2 J R J 



1 



xR 2 sm 2 e . 



(35) 



With this, the non- vanishing components of the Ricci tensor are 
(after some calculations using also the definition of the Ricci tensor in 
terms of the Christoffel symbols, identical to standard GR) 



^33 



G 



It I /2 

2 9 + 4 9 



-f + - /' 2 + JL f > 

2 J 4 J 2R J 

-f"+- f' 2 + JL/' 
2 J 4 J 2R J 

xR 2 sm 2 e . 



G 



G 



-g" + -g' 2 

2 4 

-g" + -g' 2 

2 4 



R 2 



(36) 



All other components are zero. 

The last equations were just copied from [29] . 

To obtain the tensor component we need the expression for the 
metric tensor and its inverse. We have 



g^v 



and 



1JLV 



f 1 
















-e G 
















-e G R 2 







V o 








-e G R 2 sin 2 6 


) 


/ i 
















-e~ G 














e- G 



















V o 








e~ G 


) 



(37) 



(38) 
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With this we get (K^ = g vp 7L w ) 



^ = n = (\s"+-/ 2 ) 

(39) 

The Riemann curvature is then 

n = 3(9" + 9 ,2 )-2e- G ^f" + ^ + ^A (40) 

Denoting the energy momentum tensor by T{? and exploiting the 
above results, the equations of motion are 



8lTK, 



8ttk_ 



~ 2~ J 3 



c 
8ttk 
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V 2 2R 



, \i ^ v 



A 9 



„'2 
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-off, f , 

6 iT^' 



3 /2 
I 9 



+ 6^- 



(41) 



The f'(R) refers to the derivative with respect to R, while g'(X°) 
refers to the derivative with respect to X°. The ^ functions appear 
due to the new variational principle. In [29] there appears instead the 
cosmological constant A. In principle we can add such a constant, 
too. However, one of the reasons not to do so, is that the £ functions 
will reproduce such an effect. Inspecting Eq. (14ip one can identify 
as additional diagonal contributions to the energy-momentum tensor, 
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times cr_ . Transferring £o i n the first equation of (|41|) to the left hand 
side and factorizing we can associate to £o the energy density 

PK = £-Jo*- • (42) 

Later, when we project to the pseudo-real part, this will give a con- 
tribution to the pseudo-real energy density, associated to dark energy. 
This fact will be useful in understanding the results obtained further 
below. 

4 Solving the Equations of Motion 

Homogeneity of the matter distribution requires that 



Tl = Tf = Tf . (43) 

Due to the similarity to the Tt, the same argument can be used for 
the £ functions, giving 



6 = 6 = 6 • (44) 

The £fc (k = 1, 2, 3) may be functions of the time. 

Taking an appropriate difference of the equation for T% with Tf, 
leads to the equation (subtract second equation in (j^Tj) from the third 
one in ([4"T])) 



f "~\ {f)2 ~R = ° ' (45) 

i.e, the same equation as given in [29]. The solution is also supplied 
(as proposed in [29] by direct insertion): 



1 - f R 2 



2 ' 



(46) 



with a and b as constants. Finally, the length square element takes 



the form, after some redefinitions (| ab 



14 



du 2 



(dX 



0\2 



1 



(47) 



It is exactly of the same form as in standard GR, with the difference 
that the coordinates are now pseudo-complex. This is in distinction to 
the pseudo-complex Schwarzschild metric [2D], where the differences 
appear already in the functional form of the metric. The k acquires 
the values k = 0, ±1, while Rq is a constant, related to a and b via 
| ab |= I/Rq- The function g(X°) is yet undetermined. The /c- values 
of 0, ±1 can be used to model different universes. However, as we 
know now, the k = case is the only one which is consistent with 
experiment [28] . This is the case we will finally study explicitly. It 
may serve as an example for other studies. 

We will now use the co- moving pseudo-complex coordinates [29], 
i.e., X° = 1 and X 1 = X 2 = X 3 = 0, where the dot refers to the 
derivative with respect to the eigen-time. 

The energy-momentum tensor takes the form 



(I* 



(P 



\ 



p 



\ 



/ 



(48) 



which is quite standard. The p is the matter density and p the pres- 
sure, assumed here to be pseudo-real, though, in general they can be 
pseudo-complex. We nevertheless will already restrict to pseudo-real 
functions, as they should be. 

The relevant functions in the length element take the form (see 
Eq. (13.3) of [29]) 



MX ) 
J{R) 



0\2 



R{X») 



Rl(l + kR 2 /{ARl)Y 
R{tf 

1 

^(1 + ^/(4^))' 



(49) 



which are directly obtained from the expression of the length square 
element, with some redefinitions. The R should not be confused with 
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the pseudo-complex radius variable R, nor with the Riemann curva- 
ture 1Z. R is an object with a length unit and is interpreted as the 
pseudo-complex radius of the universe. 

From now on, let us substitute X° by its pseudo-real part ct. For 
example R(X°) will be written as R(t). The derivative with respect to 
X is converted into a derivation with respect to ct, i.e., = --^ 

= BL 

c 

We part now from the above equations of motion (I41h . The ex- 
pressions in the functions / and g and their derivatives can be re 
expressed in terms of the variable R using (j49]). For example / = 
-InRl - 21n (l + kR 2 /(ARl)) and g = 2kiR. This and their deriva- 
tives have to be inserted into (14ip . Using the symmetry conditions 
of homogeneity (|43p and (|44p and the form of the energy-momentum 

tensor (|48p , the equations of motion acquire the form (remember that 
Tfl _ dR\ 

K - sr) 



8lTK 
— P 

cr 

8-KK p 

c 2 c 2 



-toe- + 



3k 3 R'{tf 



R{tY 
k 



+ 



my 



R'{tf 2R"(t) 



+ 



R(t) 2 c 2 R{t) 2 c 2 R(t) 



(50) 



The details of intermediate steps can be directly copied from any book 
on General Relativity, e.g. [29]. Instead of four equations we have only 
two, due to the symmetry conditions (|43p and (|44p . The prime refers 
now to the derivative with respect to the time t. 

Here, we will obtain one of our results, which can be extracted 
without detailed knowledge of the functions. It will be useful in or- 
der to understand the ansatz in the relation between £i and £o further 
below (see Eqs. ([541 and ([62]) ): Assuming that the density p and the 
pressure p are pseudo-real quantities, the cr_ component of Eq. ([50]) 
tells us that 



4i 



3 A; 



k 



+ 



+ 



3 R'_(t) 2 
c 2 R_(t) 2 

R'_{t) 2 2R"_{t) 



+ 



R_(t) 2 c 2 R.(t) 2 c 2 R^(t) 



(51) 



Here R-(t) is the 0"_ component of the radius of the universe. In 
order to get an idea what this implies for £o and £i, without having 
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to solve the problem, we can use the experimental result (R r is the 
pseudo-real component of the universe) 



^ = H with H' « 1 , (52) 

where the prime refers to the derivative with respect to time and H is 
the Hubble constant. Because R = R r + IRi, Ri being the pseudo- 
imaginary component of R, and I is the length parameter of the theory, 
which is extremely small (see (I)), we also can assume that R ~ R r 
and, because Ri = | (R+ — R-), we can set R± =ss R r . Using this 

R 

we can approximately write, assuming a nearly constant H = W 1 , 



R" R" R' R' R' 

= = (InRrJ^ = IHKRr)]'^ 

= [lnH + lnR,]'^ 

= H 2 . (53) 



With that, utilizing ([5T]) and k = 0, we can write the and £i 
approximately as 




e ~ 4 h2 ~ • ( 54 ) 

Further below we will see that this exactly corresponds to the case of 
a cosmological constant not changing with the redshift. This would 
be an exact result, if H is constant all over the history of the universe. 
Knowing that H is changing in time, implies that there must be a 
dependence on the radius of the universe, i.e., the redshift z. For that 
we have to solve the equation of motion exactly. As we will see further 
below, it is not easy to get the exact form of £o and £i but rather the 
use of a parametrization is appropriate. 

Let us now continue to solve the pc-RW model: 

Taking again appropriate linear combinations of (|50p (the first 
equation of (|50l) plus three times the second equation of (f50l) and 
the first equation of (|50p plus the second one) , we get new versions of 
the form 
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4-7TK ( 3p 



- (36-6) C7- 



3R" 
c 2 R 



A-7TK 



Using that 



P + 



P 



(6 - Co) a- + r2 



k R' 2 
+ 



RR" 



(55) 



RR" - R' 2 d_ (rT 

c 2 R 2 ~ dt \c 2 R / 

we arrive at the equation 



(56) 



d R' 



1 



\c 2 R) 2 VS1 su/ " ' R 2 

Differentiation of the first equation in (|50p with respect to time 



gives 



87TK dp 



d£ 6/c „, 6R' d ( 1 R' 



dt R 3 



R dt \c 2 R 



(58) 



Substituting (i57|) into (58) and multiplying the result by ^—;R 3 , yields 



R 



dp 
dt 



8ttk 



_ R 3_d^ +3R 2_ R ,_ {Ci _^ 



(J - 



3R 2 R\p + ^) . 



P 



(59) 



~~^r" — . Shifting the last term of this equation to 
the left hand side leads to 



Note that 3R 2 R' 



d / „q\ v dR 3 



8ttk 



dR J 
~dt 



(6 - &) - R 



3 *tt 

dt 



a- 



(60) 
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Identifying the mass within a given volume of the universe by M = pV, 
with V as a given volume, the last equation can be written as 



dM p dV 

~~dT + 



c 2 \dV. 




(61) 



8ttk, dt 



This is a local energy balance! In order to maintain local energy con- 
servation, we have to require that the right hand side is zero. This 
leaves us with the condition 



Any solution for £o an d £1 has to fulfill this differential equation. The 
negative index of V refers to the fact that the equation holds in the cr_ 
component. Using £1 = £o leads to ^ = 0, or £o = £1 = A = const. 
I.e., for this case we recover the model with a cosmological constant 
not changing with time. This equation is not sufficient to solve for £o 
and £1; in fact one condition is missing. 

The first equation in (|55p has the usual interpretation when £o = £1 
= 0. Then the left hand side is the sum of two positive quantities, the 
density and the pressure. The right hand side of ()55[) is proportional 
to the acceleration R" of the radius of the universe, R, multiplied by 
(-1). This equation tells us that the acceleration of R has to be nega- 
tive, i.e., we get a de-acceleration. In contrast, in the pseudo- complex 
description there is an additional term ^ (3£i — £o) °~- present, which 
might be positive. Transferring it to the left hand side may give in to- 
tal a negative function in time, i.e., depending of the functional form 
of £o and £i in time, an accelerated phase may be reproduced or not. 

Let us see whether we can get also acceleration, i.e., that 
R" > in the pseudo-complex version of GR: 

Using the left hand side of Eq. (|60j) (the right hand side is set to 
zero as argued below Eq. (|6ip ). we obtain, after multiplying with dt, 



dt 



d(lnR z _) 
dt 



(6 - £o) 



(62) 



R 3 dp + 3R 2 pdR + ^3R 2 dR 



= 



(63) 



c 



Dividing by 3i? 3 



we obtain 








(64) 



19 



Finally, dividing by (/?+$■) yields 




Now we have to make an assumption on the equation of state] 
This is a delicate part and the results can change, depending on 
which equation of state we take. The equation of state may also 
depend on different time epochs. The basic assumptions are that i) 
the distribution of the mass in the universe can be treated as an ideal 
gas, dust or radiation, the mass being equally distributed (this is only 
approximately true) . The equation of state is 



p = ap , (66) 

where p is the energy density and a is zero for a model with dust, | 
for a classical ideal gas and | for a relativistic ideal gas (radiation). 
With this, (|65p can be solved with the solution 

p = P0 R' 3{1+ ^ , (67) 

where the po is a pseudo-complex integration constant. Its dimension 
is density. 

This result is substituted into the first equation of (|55|) , solving for 
R", yields 

R" = C ^-^.-^(i + ^)poR^ . 

6 6 c z 

(68) 

We will also need the relation 



dlnV- 1 dV- Id, 3R' 

. . DO 

dt V- dt R 3 _dt ~ R 
dhiR 3 _ 



Now we remember our former result that £i has to be approxi- 
mately equal to £o (Eq. ([54])). Due to this we can assume that the 
following relation also holds approximately: 
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6 = fco 



(70) 



where (3 is an additional parameter of the theory, describing the de- 
viation from a constant Hubble parameter H. In principle, one can 
also use a power expansion of £i in terms of £oi which would only 
introduce more parameters. The f3 will later be related to observable 
quantities, like the Hubble constant and the deceleration parameter. 
Eq. ()70p gives us the missing condition, with the prize of having to 
introduce an additional parameter. Another possibility is to use the 
approximate expression of £o m terms of the ratio (R r /R r ), which 
gives £ = (3/c 2 )H, and use experimental observations for H. 
Using (|62p . we obtain for the differential equation for £o 



dt 



d(lnR 3 



dt 

d(lngg^) 
dt 



Co 



(71) 



with the solution 



£o = AJ?i (/3-1) . (72) 

This leaves us with the two, yet undetermined, parameters A and (3. 
There are many different scenarios: 

i) f3 = 1: Then £i = £o = A is constant. 

ii) (3 7^ 0: This will lead (see further below) to de-accelerated and ac- 
celerated systems, depending on the value of (3. Also the acceleration 
as a function of the radius of the universe (which can be correlated to 
time of evolution) depends on (3 and A. 

The real part of ([68]) is obtained by R'{ = \ {R'[ + R'L). Because 
the minimal length scale I is extremely small, we can assume that 
-R+ ~ R ~ R r . The real part of (|68|) is obtained by summing the 
cr+ and <r_ components and dividing the result by 2. We will also 
assume that a + ~ a„ = a, which is reasonable because the a relates 
the pressure and the density, which are both pseudo-real. Using also 
(|70p and (|72p gives the final form of the equation of motion for the 
radius of the universe 
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c 2 



R » = Lfip _ i)Ai^-D+i 

4-7tk /h , 3a x „-3(i+-%)+i 



(l + ^PoRr cI . (73) 



We have assumed that the density is real. The first term comes from 
the ^-functions. 



5 Consequences 

In this section we shall discuss the consequences of the important 
result fl73j) . 

When (3 = 1 (cosmological constant), the sign of the first term 
in (I73p is positive and contributes to the acceleration of the universe. 
The acceleration increases with the radius of the universe. For a gen- 
eral /3, the acceleration is positive, as long as (3 > g, it is negative 
(deceleration) for (3 < ^. For (3 = | no additional acceleration nor 
deceleration takes place. The last term in (|73p is always negative, i.e., 
it represents a contribution which contributes to the deceleration of 
the universe. The behavior of how the accelerating term behaves as a 
function in Rr is also determined by (3. If the exponent of R r is pos- 
itive, the acceleration increases with R r , if (3 > |, while it decreases 
with R r for (3 < |. 

The solution (|T3[) leaves space for a number of different possible 
scenarios. In order to proceed, we will make the following assumption: 
For simplicity, we assume as before that the parameter a and the 
density po are pseudo-real, i.e., a + = a_ = a and po+ = Po~ = P- in 
this case the solution simplifies to 



It; = £1(3,9 - l)AflW- 1,+1 

-^( 1 + ^)p ft- 3(1 ^ ,+ I . (74) 
3 c z 

With no dark energy, A = 0, using (fT2"j) . we have = Co = and this 
equation reduces to the one in |29] , 

Let us now discuss several particular values of f3. For that purpose 
we define 
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c 2 A 

A = • (75) 



Note that according to (|42D the A is proportional to p\, the density of 
the dark energy, with the same proportionality factor. Here po is the 
mass density of the universe. Both densities are of the same order, 
implying that A is of the order of 1. 
We can now rewrite (1741) into 



Rj, 



= A(3 i 9-l)JI^- 1 )+ 1 -(l + ^)jir 3C1+ * )+1 -(76) 

In what follows, we discuss the case of dust dominated universe, 
i.e., a = 0. For the case of a relativistic ideal gas a = |, the re- 
sults show the same characteristics. We will take arbitrarily different 
values of /3, which are chosen such that we will have the case of the 
cosmological constant A, a case which will represent the solution of the 
big-rip-off and two new solutions. These solutions are not necessarily 
represented in nature, i.e., /3 might have a different intermediate value 
as those in the examples. With this we get for 
a) = 1: (Ci = Co = A) 



R' T = 2KR r -R~ 2 . (77) 

The universe is accelerated by the first contribution and decelerated 
by the second one. For small R r the universe is decelerated. For large 
R, the first term starts to dominate and the universe is from then on 
accelerated. The turning point is at 



R r « 0.79/A^ . (78) 

If we set the radius of today at Rq = 1, a common definition of scale 
for the present epoch, the result implies that for A = 1 acceleration 
did set in after the universe passed 80 percent of its radius. This case 
corresponds to a constant cosmological function A. 
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— 3AR 2 — R r " 



(79) 



In this case, the acceleration increases with the second power in R r , 
stronger than only with a cosmological constant. The £o function (|72p 
is then given by AR-, i.e., the dark energy density, represented by £o> 
increases with the radius of the universe.. This is like the big rip-off, 
which is discussed in the literature. The break-even point, i.e. when 
acceleration is equal to deceleration, is reached for 



Rr = 1/(3A)3 . (80) 

For A = 1 the break-even point is reached when the universe is 1/3 of 
its present radius, thus, earlier than in case a). 

c) /3 = 2' Remember that a = (dust dominated universe)! Then, 
from ([76]) we get 



R" = ^AR r * - R~ 2 . (81) 

In this situation, the dark energy behaves as (use Bq. <E>) Po = jfr, 

i.e., the density of the dark energy decreases with the radius (time) of 
the universe. 

This is really a new solution! The accelerating and the decelerating 
parts are decreasing with the size of the universe, but at a different 
rate. For small R^, the second term dominates and the universe is 
decelerated, while for sufficient large Rr the first, accelerating, term 
dominates and the universe is accelerated! The break-even point is at 



R r « 2§/A3 , (82) 

i.e., for A = 1 the universe at this break-even point will be at about 
25 « 1.59 times of its present radius. Of course, this can be changed 
using different values of A. For A = 3 the break-even point is at 76 
percent of the radius of the universe. This case is plotted in figure 1. 
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The universe starts to be accelerated after having reached R r = 0.76 
(units in Rro). After that the acceleration increases. However, having 
reached the radius R r ~ 1.9, i.e., nearly twice the actual radius of 
the universe, it reaches a maximum and after that the acceleration 
is decreasing, reaching asymptotically zero. This universe will never 
collapse but reach an asymptotically non-accelerating state 



4 - 



2 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
radius 



Figure 1: Dependence of the scaled acceleration as a function of R r , the 
radius of the universe, for (3 = \. In this figure A = 3. The maximum of 

2 

this function can be deduced from f lHTj) . giving i? max = (s/A^ 3 = 1-923. The 
maximum can be barely seen in the figure due to the extreme slow decrease 
of the function. The question is also: Where are we now? Before or after the 
maximum? 

d) j5 = |: Remember that a = (dust dominated universe)! Then, 
from ([75]) we get 



R r — A — R r ^ . (83) 
This is also a new solution. The break-even point is now at 

R r w l/y/X . (84) 
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This solution is also special in the sense that the asymptotic accel- 
eration of the universe is constant (R" = A). Using (f72|) leads to 
the dependence po = A/R r of the dark energy on the radius of the 
universe, i.e., it also decreases with the radius (time) of the universe. 



4 - 
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Figure 2: Dependence of the scaled acceleration as a function in R r , 
radius of the universe, for (3 = |. In this figure, the A = 4. 

In all cases A can be fitted to the observation at which stage the 
net acceleration did set in, overcoming the deceleration term in (|74D . 
The new part here is that other solutions exist than the standard ones: 

i) There is the possibility of a constant asymptotic acceleration. 

ii) In another solution, the expansion of the universe, after its decel- 
erating period, gets accelerated. The accelerations reach a maximum 
and vanish asymptotically. In this case the universe approaches, for 
large times, an ever expanding, non-accelerating phase. 

iii) In all cases, the universe is first decelerated and after a so-called 
break-even point it starts to accelerate. 

iv) Of course, all standard solutions are obtained (cosmological con- 
stant and rip-off). 

In order to calculate numerically observable consequences, we have 
to know the exact form of the functions, which we were unable to 
deduce from first principles. One possibility is to use the calculated 
distribution of dark energy, as for example done in [22]. An alterna- 
tive is to use the parametrization given in (|70p . This implies the use 
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of an additional parameter (/3) and is equivalent to known consider- 
ations in the literature [28]. The acceleration in each solution is a 
consequence of the functions. As discussed above, they represent 
contributions to the energy-momentum tensor, equivalent to 
the dark energy. In the model considered, this dark energy 
is in general not a constant but may vary in time, i.e., with 
the radius of the universe. 



Extraction of j3: 



We can try to connect the value of f3 to observable quantities. For 
that, we start from Eq. (|53|) . without the approximation in the last 
line. We get 



R 1 '. 
R r 



H' R' r 
H R r 



El 



+ H 



R r 

H = H' + H 2 



(85) 



Substituting this into the expression for £1 (Eq. ([5"T]) ). setting k = 0, 
we obtain 



1 H 2 + ^(H' + H' 



—^H 2 H — kH' 



= fco , (86) 

where we have used that £i = (3£q. Using = ^H 2 (see (|MD ) and 
solving for (3, we obtain the final result 



" - 1 + ■ < 87 » 

We obtained further above that for (3 > g the universe is accel- 
erated after a given radius. This corresponds to M% > — 1. In order 
to get deeper insight, we use the deceleration parameter, which is a 
measure whether the universe is accelerated or decelerated, depending 
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on the sign of this parameter. The deceleration parameter is defined 

The universe is 

1 „„ H> 



as [29] q = — ^^M T , which gives q 



accelerated when (3 < |, or ^ < — 1, or q > 0. The acceleration 
increases with R r when /9>|, or -j^ > — |, or q 1 < — |. 

In conclusion, a measurement of the change of the Hubble constant 
with time will lead to a determination of the parameter /3 as a function 
of time. Though, the last considerations clarify the role of (3, we are 
suffering still by the problem that we have to know the solution of 
R' 

H = Wj. This can be done up to now only through the experimental 
measurement of the Hubble parameter H. 

A model including dust and radiation, k=0: 

Up to now, we did only consider one density component (dust 
or radiation) and the pseudo-complex contribution, given by the 
functions. Realistic models involve both components, as can be seen 
in |28j . Expressing the ratio of the radii Rj-q and -R r , the present 
radius of the universe and the one at a redshift z, respectively, in 
terms of the redshift, we obtain [291 EH] 



^ = (1 + *) , (88) 
it, 

where the index r refers to the real value of the radius. We obtain for 
the ratio of the velocity and the the radius of the universe [28] 



^ = h = H${n d (i + z) 3 + n r (i + z) 4 + n A f( z )} 

(89) 

where the index d refers to the dust part and the index r to the 
radiation part. We do not include the contribution due to k ^ 0, 
because we consider a flat universe. The factor Hq is the square of 
the present Hubble constant. Using our previous result, the function 
f(z) is given by 



f(z) = (1 + z) 3 ^" 1 ) = (l + zf( 1+w *> , (90) 

(/3 = 2 + w) where we made a connection to the notation used in [28]. 
Our result states that when the Hubble constant changes in time, 
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there must be a deviation from f3 = 1, which corresponds to to = —1, 
the case of a cosmological constant, constant in time. The deviation 
from (3 = 1 cannot be large when H', the time derivative of the Hubble 
constant, is small. Up to now, we only find a parametrization of the 
functions in terms of the Hubble parameter or the parameter /3. The 
deeper origin of the value of can probably explained by fundamental 
theories like string theories. 

Note, that the relation (3 = 2 + w with (|87j) gives a relation of w to 
the Hubble parameter and its derivative in time, which makes definite 
predictions on w, once H and H' are known. To our knowledge, this 
is not presented elsewhere. 

6 Conclusions 

We have applied the pseudo-complex formalism to extend the Robertson- 
Walker model of the universe to the pc-RW model. The main results 
are that 

1) The model introduces automatically a contribution which is equal 
to the cosmological constant or dark energy which may depend on the 
radius of the universe. 

2) The cosmological "constant" is a constant when the Hubble con- 
stant is constant too. When the Hubble constant changes slightly with 
time, our model predicts deviations from the cosmological constant, 
depending on the redshift (time of expansion). The amount of devia- 
tions depends on the exact form of the functions. 

3) The deviation can be obtained, once the radius of the universe, as 
a function of time, is known. Within our theory, we obtain several 
possible dependencies of the dark energy as a function in the radius 
of the universe, depending on the parameter j3. 

4) We also obtained several possible evolutions of the universe. Be- 
sides the solution of a constant dark energy density and the rip-off 
scenario, we also obtained solutions where the acceleration tends for 
infinite time towards zero or a constant value (see Figs. [Hand [2]). 

5) We obtained a relation between w = f3 — 2 and H and H'. Once H 
and H' are known, the w value can be deduced. 

The origin of the functions might have a deeper microscopic 
origin, which we do not explore here. Probably, only such a deeper 
microscopic understanding will fix the dependence of on the ra- 
dius of the universe (see for example [22]). Nevertheless, the classical 
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picture presented here enlightens and simplifies the description of dif- 
ferent possible evolution scenarios of the universe. 

We have not yet investigated the role of the minimal length param- 
eter I, which also appears in the pseudo-complex formulation. In field 
theory its function is to render the theory regularized [30] . We suspect 
that this also happens in the pseudo-complex formulation of General 
Relativity and might give a hint on how to quantize this theory. In 
a future publication we intend to investigate the role of the minimal 
length scale I. 

We saw that the modified variational principle 5S e P° has impor- 
tant consequences as the appearance of dark energy. It also provides 
a simpler description of effects of the dark energy, obtained via quite 
involved numerical calculations, as for example in |22j . These features 
are a hint that the variational principle has to be probably modified 
as proposed. 
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